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' Abstract. We prove a Wegner estimate for a large class of multi- 

I , particle Anderson Hamiltonians on the lattice. These estimates will 

• allow us to prove Anderson localization for such systems. A detailed 

| proof of localization will be given in a subsequent paper. 



' 1. Introduction 

^2 \ Wegner estimates originate in the famous paper [10] . There, Wegner 

proved among other things that the integrated density of states for the An- 
derson Hamiltonian has a bounded density provided the probability distri- 
| bution of the random potential itself has a bounded density. This implies in 

particular an upper bound on the probability that an Anderson Hamiltonian 
on a finite box has eigenvalues close to a given energy E. 

Wegner's estimate play a key role in the Multiscale method to prove 
Anderson localization (see e. g. [5] or [4]). Only recently Bourgain and 
rN \ Kenig [TJ proved Anderson localization for a Bernoulli model without an 

a priori Wegner estimate; they prove a Wegner-type estimate inductively 
within the Multiscale scheme. 



Wegner's original work was restricted to lattice models. However, the es- 
timate was also proven for the continuum (see [3] for a recent rather optimal 
result and [9] for a review on this subject). 

In this note we prove a Wegner estimate for a multi-particle Anderson 
model. In a subsequent paper we will also do multiscale analysis for this 
model. The first Wegner estimate for a multi-particle random Hamiltonian 
was proved by Zenk [11] . Chulaevsky and Suhov [2] develop a multiscale 
analysis for certain (1-d) two body Hamiltonians. The Wegner estimate 
in this paper requires strong conditions on the probability density of the 
random potential (e. g. analyticity) . It was one of the motivations of the 
present note to avoid these strong assumptions. 
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The method of proof applied here is close to Wegner's original idea and 
was developed from the paper [6]. Note that there is a refinement of this 
method by Stollmann [8] which is likely to work in the multi-particle case 
as well. 

We note that the method presented in this paper will also work for alloy- 
type models in the continuous case. The necessary changes can be read off 
from the paper [6]. However, in the continuous case we get the volume factor 
of the bound with an exponent 2. This suffices to do a multiscale analysis, 
but it gives no result for the regularity of the integrated density of states. 

2. Models and Results 

We will deal with a system of N interacting particles on a lattice Z d . We 
consider these particles on the full Hilbert space, disregarding Fermionic or 
Bosonic symmetry. Physically speaking we deal with distinguishable par- 
ticle. Since the full Hilbert space is a direct sum of the irreducible sub- 
spaces with respect to iSjv-symmetry (including the totally symmetric and 
the totally antisymmetric subspaces) the Wegner estimates for Fermions and 
Bosons follow immediately from the result on the full space. 

The one-particle Hilbert space we consider is £ 2 (Z d ) and the Hilbert 
space for N (distinguishable) particles is consequently i 2 {'L Nd ). Any (bounded) 
operator A on these Hilbert spaces is uniquely defined through its matrix 
elements A(x,y) = (5 x ,A5 y ) where 5 Z is the vector in £ 2 with component 
1 at lattice site z and otherwise. 

We write the lattice site x £ Z Nd as x = (x±, . . . , xjv), where Xi G Z d 
denotes the the coordinates of the i th particle. 

Each particle (with coordinates £) is subject to a random potential f w (^) 
which is the same for all particles. The random potential v^^) consists 
of independent identically distributed random variables. Throughout we 
assume that the distribution of the v(£) has a bounded density p{v). We 
denote the underlying probability measure by P and the expectation with 
respect to P by E. 

The kinetic energy operator for one particle is given by: 

hou(0 = u (Z + n ) (2.1) 

H=l, n&L d 

the single particle Hamiltonian is consequently: 

K = h Q + v u (2.2) 

If h is a one-particle operator acting in £ 2 {7 l d ) we denote by the 

corresponding operator on £ 2 (Z Nd ) acting on the i th particle only, more 
precisely: If h has matrix elements /i(£, rj) then: 

/i w u(xi, . . . ,x n ) = ^2 h(xi,rf)u(xi,... ,Xi-i,r},x i+ i,...,x N ) (2.3) 
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In other words, 

= 1P(ZJ-) ® • • • ® lp(Zd) ® & ® ^(Z*) ® • • • ® lf2 (Z d) (2.4) 
V v ' V v ' 

i—1 times n—i—l times 

The N-particle Hamiltonian without interaction is defined by: 

N 

H^o = Yl hil) ( 2 - 5 ) 
i=i 

The interaction term U can be a rather general function on Z Nd . We as- 
sume it to be bounded for simplicity. We also suppose that U is a determin- 
istic function, it would be sufficient for our purpose to have U independent 
of Vw In most cases U is a pair potential of the form U(x) = u(xi-Xj). 

The N-particle Hamiltonian with interaction Uis then given by: 

H u ,u = Hu, + U (2.6) 
We will deal with this operator restricted to a bounded (hence finite) 
domain A. The number of elements of A will be denoted by | A |. 
We call a subset R of 7L d a rectangle if 

R = g Z d | L v < i v < M v for v = 1 . . . N } (2.7) 

A rectangular domain in r L Nd is a set A of the form: 

A = Ai x A 2 x...xAjv (2.8) 

where the Aj are rectangles in Z d . We use the notation IT (A) = A». We 
call a rectangular domain A regular if for all i, j = 1, . . . N either Aj n Aj = 
or Aj = Aj. 

For any subset A of Ij Nd we define the operator = v by its matrix 
elements: 

H A (x,y) = H^u{x,y) forx,yeA (2.9) 

The main result of this note is the following Wegner estimate for multi- 
particle operators: 

Theorem 2.1. If A is a regular rectangular domain then 

P ( dist ( (7 (F A ) , E) < ft) < C || p ||oo | A| k (2.10) 

The assumption of regularity of the set A can be avoided. However, 
the proof is more transparent with this assumption. The proof of Anderson 
localization by multiscale analysis which we will present in a forthcoming 
paper will deal with regular domains only. 

The proof of Theorem 12.11 implies also that the integrated density of 
states has a bounded density. This result can also read off from the explicitly 
known form of the integrated density of states (see [7]). 
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3. Proof 

We prove Theorem 12.11 Let A = Ai x A2 • • • x A^v- We may assume that: 

Ai = A 2 = • • • = A K (3.1) 
and Ai n Aj = for all i > K (3.2) 

We denote the eigenvalues of H A by E n = E n {H A ). We order them so 
that E\ < E2 < . . . and repeat any eigenvalue according to its multiplicity. 
The eigenvalue counting function is denoted by: 

N(H A ,E) = #{E n (H A )<E} (3.3) 

We will need the following Lemma: 

Lemma 3.1. Suppose \3. 1\) and \3. 2ty hold. Denote by v(£) the value of 
the random potential v u evaluated at the lattice site £ € Z, d . Then: 

Proof. Set V(x) = J2?=i vfe) Then for ^gAj: 

J^(xi,...,x n ) = k*i ( 3 - 5 ) 
Hence for each [x\, . . . , xjy) £ A we have 

dv K 
S ^7?T (:El '---' X7v) = ^ S s t** = K (3 - 6) 

IgAi £eAi i=i 

Let us denote by ^ n the normalized eigenfunction of H A for the eigen- 
value E n = E n (H A ). The Feynman-Hellman- Theorem tells us that: 

dEn > , dv 1 \ sr 1 / < \ 12 aF /o 7 , 



Thus from (13.61) we obtain: 



E^-Ei*.wi'(E^) <-» 

= K ^ \ M*)\ 2 =K (3.9) 
zgA 

since ip n is normalized. □ 

Let (p be an increasing C°° — function on R, < ip < 1 with 93 = 1 on 
(k, 00) and </? = on (—00, —n). 



A MULTI-PARTICLE WEGNER ESTIMATE 

Then: 



dist(a(# A ) , E) <k) (3.10) 

< E ^N(H a ,E + k) - N(H A ,E- «)) (3.11) 
= E(tr ( X{E _ KtE+K] (H A ))) (3.12) 

< E ( tr ( ^(i/ A -E + 2k) -ip (H A - £- 2kJJ (3.13) 

<E^ tr (V(# A - £ +t)) cft^ (3.14) 
by Lemma 13. II 

4 E £ « (^»< ffA > " * +') E ^) * (3.15) 

.a, ,df(E n (H^) -E +t) s 

^E /. E ■( v a„ K) < 3 - 16 > 



2K ^GAj 



Since E is a product measure we can split it into an integration over v(£) 
which we write as J ■ p(v) dv and the expectation with respect to the other 
random variables, which expectation we denote as . 
With this notation (|3.16p equals: 

f2« - r t? I trA 



E / " * E K m ( / MEm Z a E + t) MO ) (3.17) 

E/>'E^(/ a ^r +t) ^) <3 - i8> 

n J—AK tel. ^ V>/ 



€eAx 

(•2k 

By the fundamental theorem of calculus we have: 



dip (E n (H A ) - E + t) , , , 

= v(E»(fl„ A ,0=„«) - E + 1) - f - £ + t) (3.20) 

where iJ (resp. il,,-, ■ ) denotes the operator H A with the po- 

tential v(£) set to its maximal (resp. minimal) value, i. e. with = 
sup (supp(p)) or v(£) = inf (supp(p)). Note that we include the cases 
sup (supp(p)) = oo and inf (supp(/?)) = — oo. 

Changing i>(£) from its minimal to its maximal value is a (positive) 
perturbation of rank at most M = Kjxjj- Thus: 

E n{H A ^- )=min ) < E n {H A ^ =max ) < E n+M (H A ^ =rnin ) (3.21) 
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To estimate (I3.18P we use the following simple Lemma: 

Lemma 3.2. Let if be a non decreasing function on M. with < ip < 1. 
// a n and b n are non decreasing sequences satisfying a n < b n < a n+ M for all 
n, then: 

( (p(bn) - <p(a n ) ) < M (3.22) 

n 

Combining the above estimates we get: 

dist ( a (H A ) ,E) <k) (3.23) 



CeAi 



< IIP II ^ dtJ2 E- {0 Y, (v{E n {H^ )=max ) - ip{E n {H^ )=mm )) 

J -2k ^ eAi n 

< \\p\\ Ak |A| (3.25) 
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